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Abstract: We study a coupled dark energy-dark matter model in which the energy- 
K> , momentum exchange is proportional to the Hubble expansion rate. The inclusion of 

its perturbation is required by gauge invariance. We derive the linear perturbation 
equations for the gauge invariant energy density contrast and velocity of the coupled 
fluids, and we determine the initial conditions. The latter turn out to be adiabatic for 
dark energy, when assuming adiabatic initial conditions for all the standard fluids. 
We perform a full Monte Carlo Markov Chain likelihood analysis of the model, using 
WMAP 7-year data. 



1. Introduction 

The true substance of dark energy and dark matter is unknown although it should 
account for about 95% of the matter-energy content of our universe today [1]. While 
the couplings of dark fluids to photons and normal matter are severely constrained [2] , 
nothing prevents dark matter-dark energy interactions [3-21]. At the level of the 
background evolution equations, it is customary to parametrize the coupling between 
the two dark sectors [22] as: 

Pdm + ^y.Pdm = O.Qdm ' (1-1) 

Pde + ^T-ipdei^ +w) = aQ^^ , (1.2) 

where pdm, Pde denote the dark matter and dark energy energy densities, respectively, 
and Q^^ = —Q^e encodes the coupling between those two dark sectors and drives 
the energy exchange between them. The dot indicates derivative with respect to the 
conformal time dr = dt/a, with "H = d/a = aH denoting the background expansion 
rate, while w = Wde = Pde/ Pde stands for the background dark energy equation of 
state and pressureless dark matter is assumed: Wdm = Pdm/ Pdm = 0. From now on, 
barred quantities are to be considered as the background quantities. 

The initial conditions for the several components populating the early universe 
have been explored to a large extent. They were first analyzed for all cosmic flu- 
ids but dark energy (see e.g. Ref. [23] and references therein), with the result that 
adiabatic initial conditions were one possibility. It was also noticed that the choice 
of gauge could be a delicate issue: a safe alternative proposed was to use a gauge 
invariant formalism [22,24-26]. The initial conditions for the case of dynamical dark 
energy as an uncoupled quintessence field have been also derived [27-33], including 
a gauge invariant treatment [33]: they turned out to be adiabatic if those for the 
traditional fluids were adiabatic. Furthermore, the formalism in Ref. [33] has been 
recently applied to the case of a coupled dark energy-dark matter systems which 
mimic uncoupled models at early times, both at the background and perturbation 
levels [18] for the viable parameter space: as expected, adiabatic initial conditions 
for dark energy naturally resulted then. Here we consider a different class of dark 
couplings, not negligible at early times. It is also illustrated that the gauge invariant 
formalism is particularly illuminating for the determination of the correct perturba- 
tion equations, for a general coupled theory. 

The structure of the paper is as follows. In Section ||, the notation is set and 
the gauge invariant equations -at linear order in perturbation theory- for a coupled 
fluid are derived. In particular, we study in Section ^]2| the case of a (covariant) 
dark matter-dark energy interaction proportional to the Hubble rate. In Section ^ 
following the method proposed in Ref. [33], we derive the corresponding initial con- 
ditions for dark energy. Then in Section ^ we constrain the type of coupled models 
analyzed, using several data sets. Section ^ contains the conclusions. 



- 1- 



2. Gauge invariant perturbation equations 

Following Ref. [22], the FRW metric, up to first order in perturbation theory, can be 
written as: 

g^^dx^'dx" = a^[-{l + 2A)dT^ - Bidrdx' + (7^- + 2Hij)dx'dx^] , (2.1) 

where jij is the 3D flat metric with positive signature. The perturbations A, Bi and 
Hij are functions of time and space and are in general gauge- dependent, i.e. not 
invariant under an infinitesimal coordinate transformation: 

(a;°, x') -^ (x°, x') = (x° -T,x'- U) . (2.2) 

Particularizing to the case of scalar metric perturbations, two gauge invariant quan- 
tities^ can be defined [24], the most popular being the so-called Bardeen potentials 
$B and '^B- 

In describing the evolution of a given fluid "a" , other gauge dependent quantities 
are introduced, such as the perturbed 4~-velocity and the energy-momentum tensor, 
which can be expressed by: 

< = ^(1-A<), (2.3) 

T^'' =-pa{l + 5a)K< + r!i\ (2.4) 

where f ^ is the peculiar velocity perturbation of the fluid, 5a the density perturbation 
and T^'^ the stress tensor, whose components in first order perturbation theory read 

T,° = 0, rJ, = Pavl, ^5=p4(l + ^^)^ + (7rJ)5] . (2.5) 

In what follows, we deal with the Fourier transformations of the scalar part of the 
metric and fluid perturbations. See Appendix ^ for details. In the equations above, 
vrf' and vrj denote the isotropic and anisotropic scalar pressure perturbations, respec- 
tively, while Va is the scalar part of the peculiar velocity. Associated gauge invariant 
quantities can be defined, paralleling the two gauge invariant variables for scalar 
metric perturbations. Following the notation in Ref. [33], a possible gauge invariant 
formulation for (5^, Va and the stress-tensor components vr^ and vrj is: 

(2.6) 

(2.7) 

The coefficient c^„ entering in the entropy perturbation F^ is the adiabatic sound 
speed of the fluid c\^ = pa/ Pa and Wa is the equation of state of the fluid. 

We focus next on the derivation of the gauge invariant equations for the matter 
density contrast A^ and the fluid velocity Va, for a generic coupled fluid. 
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^Thc transformation properties of the metric perturbations defined in Eq. ( |2.lD and the exphcit 
definition of the Bardeen potentials is reminded in Appendix ^ 
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2.1 Coupled fluids in general 

Consider the full (background plus perturbations) continuity equation for fluid "a" 

^,Tli'' = Q: , E^^ = 0' (2.8 



where T^'^ denotes the corresponding energy-momentum tensor and the vector Q'^ 
governs the energy-momentum transfer. The constraint on the right accounts for 
total energy-momentum conservation. Following Ref. [22], Q'^ can be written as: 



SQa 
Qa 



Qa = Q J 1 + 



with 



= Qail+ea) 



3>l = 



(2.9) 
(2.10) 



where j^ and Ea are perturbation parameters. In particular, the background contri- 
butions reduce to the coupled dark energy-dark matter case in Eqs. ( |1.1| ) and ( |1.2| ), 
for Ql^g = —Qdm- Defining for simphcity j* = pafa/a, the total coupling reads 



Q'a 



1 /^ 



■jl-(A-Ea)],QX + Para) 



(2.11) 



Let's denote by /„ the Fourier transform of the scalar part of /^. One can show that 
fa is invariant under gauge transformations, while £„ transforms as 



? =e - ^T 

-a ^a 7n 



(2.12) 



where the "hat" denotes gauge transformed quantities. This suggest a possible choice 
of gauge invariant variables for the coupling perturbation parameters, given by 



E. 



Sa- 






^ a J a ■ 



(2.13) 
(2.14) 



The gauge invariant choice in Eq. (|2.13|) is analogous to that for A^ in Eq. ( p.6|) . 

With the help of these variables, the scalar perturbation equations for the matter 
density contrast A^ and the peculiar velocity \4, for a generic coupled fluid, read: 

A, = -m [(4, - wa) A, + ^,r„] - k{l + Wa)Va + 3Hg„ [A + E^- AJ , (2.15) 



K = -H{l- 3cL) Va + - 



+ Wa 



C^^Aa + Wa\Ta- -Tia 



-snqa 



Va - k^) + -^^ 
l + WaV" 'H J ^+Wa 



^Aa 



k {^n - 3cL$s) 
(2.16) 
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where ^ is a metric gauge invariant quantity [22], whose expression is given in 
Eq. (|A.19|) . The quantity qa accounts for the energy transfer Q^ in Eqs. ( p.l5| ) 



and ( p.l6p , rescaled as follows 



^'-wk- <^-"* 



For vanishing qa and Fa, Eqs. (|2.15|) and ( p.l6| ) reduce to those in Ref. [33]. 



2.2 Coupling proportional to H 

Coupled models with a dark matter-dark energy coupling proportional to the Hubble 
expansion rate have been studied at the level of linear perturbations in several recent 
works, see for example Refs. [14-17,34]. Perturbations in the expansion rate were 
neglected, though. To analyze the issue, the results of the previous section will be 
particularized to the following coupling: 

Q:^^ = ei^p,,nL = -Q:^e- (2.18) 

Here Q'^^ is chosen parallel to the dark matter four velocity u|^^ to avoid momentum 
transfer in the dark matter rest frame [14]. The evolution equation for the dark 
matter velocity remains then equal to that of baryons, avoiding the violation of the 
weak equivalence principle. Moreover, the authors of Ref. [14] pointed out that such 
a coupled models could suffer from non-adiabatic instabilities if the coupling Qdm is 
chosen proportional to the dark matter energy density. In Ref. [15-17], it was shown 
though that such instabilities could be avoided in a minimal way choosing a coupling 
Qdm proportional to the dark energy density^. 

It is important to notice that, in order to deal with a consistent model, H in 
Eq. (|2.18| ) must denote the total expansion rate (background plus perturbations). 



H = H + 6H, while in all previous studies only the background quantity was consid- 
ered. The inclusion of SH is mandatory to preserve gauge invariance, as we proceed 
to illustrate. For the model in Eq. (p.l8|) one obtains: 

Q,„, = (H-pde, (2.19) 

edm=^-^ + ^-^^IC + 6de. (2.20) 

H Pde 

The K, term (see Eq. (|A.9|) ), represents the expansion rate perturbation, overlooked 
in all the references mentioned above. Indeed, K, depends on the time slicing, so that 
the coupling perturbation Sa gauge transforms as: 

e,-e,^5iT = |r + ^r = (Z-lc) + {5de-6de)- (2.21) 



^Would an interaction proportional to the dark matter density be studied instead, it would be 
necessary to consider a time dependent dark energy equation of state, in order to avoid early time 
instabilities, thus introducing at least one extra free parameter, see Ref. [18]. 



-4- 



To our knowledge this result was not explicitly discussed elsewhere. We will see 
in Sec. ^ that the extra contribution resulting from 5H has little quantitative im- 
pact on the physical constraints obtained from data, while being essential for gauge 
invariance. 

Before proceeding further let us comment on the covariance of the coupling of 
Eq. ( 12.181) . First of all, the dark energy density can be rewritten as pde = Tj^e'^'lf'^u^- 
Moreover, we can express the Hubble expansion rate in terms of the covariant deriva- 
tive of the four velocity defined in Eq. (|2.3| ). Indeed, it is straightforward to verify 
that the background quantity associated to m^.^ is directly proportional to the ex- 
pansion rate H. Following [22] one has: 



e, = <,. = 3H{1 + ICa) . (2.22) 



Under gauge transformations, the perturbation /C^ (associated to the a-fluid) trans- 
forms like: 



]Ca-JCa = =T (2.23) 



which is exactly what is needed to preserve the gauge invariance of the coupled model 
under study, see Eq. ( p.21| ). In the following, we will use for definiteness the total 



matter expansion rate 9^ = ui^ = 3H{1 + /C), denoting with /C the perturbation 
associated to the total fluid. Finally the coupling of Eq. ( p.l8 ) can be written in a 
covariant way as: 

Q^m = ^^ Tl^<'uf uL = -QL ■ (2.24) 



We can now particularize Eqs. ( p.l5| ) and (|2.16|) to our coupling. Expressing /C 



in terms of gauge invariant quantities one obtains: 



Ea = Ad,+ [j-^{l + WT)]VT + 2$B (2.25) 



where wt and Vr is the equation of state and velocity of the total fluid. The density 
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and velocity perturbation equations then read: 



A 



dm 



n 

•^ dm 
Ade 



-X Vdm + ^- 



Pde 



Pdv 



X 



{Ade - Adm) + Y^T 



2L 



K 



dm 



<^B + fi^n, 



(2.26) 
(2.27) 



n 



-3(4 - w)Ade -{1+w) X" Vde + 9 (1 + W) (4 - 4) Ub ' Vde 



< 



X 



Vt - 3 (4 - 4) Ub - Vde 



f2.28) 



v^ 



de 



n 



2L 



+ 



e 



l + w 



34 ]Vde-{i + 34) $B - n,u, + -^Ade + 

I l + w 

(1 + 4) Vde - Vdm - cl<^B , 



(2.29) 



the rescaled quantities V = V/x and 11 = U/x"^ were used, with x = k/Ti. In deriving 
these equations, the dark energy entropy perturbation Tde has been rewritten in terms 
of Ade, Vde and $b, see Eq. ( |A.26| ). 4 ^^^1 4 are the dark energy adiabatic sound 
speed and the rest frame sound speed, respectively. In the following we work in the 
framework of constant w, c\ = w and 4 = 1- 

3. Initial conditions 

In Ref. [33], whose gauge invariant formalism we follow, the solution of the system 
of differential equations for the perturbations is reduced to that of a simple eigenval- 
ues/eigenvectors problem: 



U' 



dU 
dlnx 



A{x)U. 



(3.1) 



Here A[x) encodes the evolution equations for all the universe components, and 

U = {Adm,Vdm,A^,V^,Ah,A^,V^,Il^,Ade,Vde} (3.2) 

is an array of gauge invariant perturbations, where the subscripts 7, b and u stand for 
photons, baryons and neutrinos, respectively. No anisotropic stress for dark energy 
and negligible anisotropic stress for photons (due to large Thompson damping) are 
assumed. 

The evolution equations for baryons, photons, and neutrinos are unaltered by 
the presence of the dark coupling and we obviate them below. In contrast, the 
dark matter and dark energy perturbation equations for the case under study are 
significantly modified. The exact form of the correspondent A{x) matrix can be 
easily derived from Eqs. ( ^.26[ )-( ^.21J ). 
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To obtain the initial conditions for cosmological perturbations, it is necessary to 
study the evolution of the several cosmic components at a very early stage, when 
the universe was radiation dominated and l-L = 1/r. One is interested in the time 
dependence of all perturbations on super-horizon scales, i.eioi x = kr <^ 1. 

3.1 The Aq matrix 

At early times x ^ 1, the A{x) matrix can be approached by a constant matrix 
Aq, if no divergence appears when taking the hma;_j.o A{x). The assumption that the 
universe is radiation dominated at early times implies wt = 1/3, pr = prad and 

n, = Pu/Prad = Ru, n^ = l-Ru and ^ = l^ oc x'^^'^+^K (3.3) 

^ ^dm Pdm 

w < —1/3 is assumed as well, in order to obtain cosmic acceleration, which implies 

that {3w + ^) can be always taken negative for ^ < 0. 

Using Eqs. (|2.26| )- (|2.29|) and taking the a; — )■ limit, the following entries in the 
Ao matrix associated to A^g and Vde result: 

'O ^(« + /?0 R-,ia + m %(a + /3e) R,Aa + m -/3 -{a + PO' 

-Cr -Ry (1 + er/4) -4i?^ (1 + er/4) - R,. (1 + Cr / 'i) -4i?. (1 + er/4) - R. J^ 1 + 2^,. 

where a = 9(1 — w^), (3 = 3(1 — w), ^^ = ^/(l + ui) and R^ = 1 — Ry. The other 
lines in the Aq matrix remain equal to the uncoupled case ones. Indeed the extra 
term in the ^dm equation proportional to ipdelpdm can be safely neglected in the 
X — ;■ approximation. We thus recover the standard non interacting dark matter 
perturbation equation in the early universe. Notice that this was also the case of the 
viable coupled model discussed in Ref. [18]. 

3.2 Adiabatic initial conditions 

Let Ui be an eigenvector of Aq with eigenvalue Aj. The solution to the system in 
Eq. (|3.1| ) can then be expressed as a linear combination of x^^Ui terms. Those cor- 
responding to the largest eigenvalues will dominate the time evolution. We checked 
that for the model under study, Eq. ( |2.24| ), the dominant modes are associated to 
Aj = values and they suffice to specify the initial conditions. The subdominant 
modes decay in time as they correspond to negative real eigenvalues^. The dominant 
eigenvalue of the evolution matrix, Aj = 0, is fourfold degenerate (as was the case for 
a universe without dynamical dark energy) and the corresponding four eigenvectors 
serve as a convenient basis to specify the initial conditions. 



■^Also, see Ref. [33] for the case of quintessence and Ref. [18] for coupled dark sectors with a 
different coupling. 
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Let us assume adiabatic initial conditions for all species but dark energy, as 
strongly constrained by WMAP data [1,35]. For each pair of components oi and 02, 
the relative entropy perturbation, S'^^qj, vanishes: 

AO AO 

Sa,a, = -^^ - ^^ = . (3.4) 

Pai/ Pai Pa2/Pa2 

For baryons, neutrinos, photons and dark matter this implies: 

(3.5) 





dm, 


= A° = -A° = 

fe 47 




from which one obtains: 










= vL 


= --PAO 


and 



no = -pa; , (3.6) 

with V = 1/(15 + 4:R^). Those are the standard adiabatic initial conditions for 
velocity perturbations and anisotropic stress. Solving the eigenvalue problem for our 
Aq matrix, it follows that dark energy also obeys adiabatic initial conditions given 

by 



i(— f 



Al='-{l + w+'-)Al (3.7) 



K.°=-^A?- (3.8) 

Consequently, adiabatic initial conditions for the matter and radiation components 
automatically imply adiabatic initial conditions for dark energy, alike to the case 
for tracking scalar quintessence [33] or those obtained for dark energy-dark matter 
couplings which do not depend explicitly on the Hubble rate^. 

As a final comment, notice that the previous results do not depend on the fact 
that we are using the expansion of the total fluid Ot (and its perturbation /C) to 
define the dark coupling in Eq. ( |2.24| ). In fact one could have used the expansion 
rate of any single specie, 6a. In that case, Eq. (|2.25|) should be replaced by: 



E, = A,e + y^« - ^(1 + '^t)Vt + 2$B (3.9) 

which implies that all the contributions of the dark coupling going as x'^Vr in 
Eqs. ( |2.26|) and ( |2.28|) should be replaced with x^Va- This does not modify the 



expression of the matrix Aq of Sec. |3Tl| , that encodes the evolution equations at early 
times, as in the limit a; —> all the x^-terms can be neglected. As a consequence our 
results of Eqs. ( |3.7| ), and ( p.8| ) would not be affected, and our conclusion on adiabatic 
initial conditions remains unchanged. 




Figure 1: Left (right) panel: la and 2a marginalized contours in the ^-Qdm.h'^ (C^w) 
plane. The contours show the current constraints from WMAP7, HST, SN, H[z) and LSS 
data taking into account the expansion rate perturbation /C. 



4. Data constraints 

In this section we briefly revisit the constraints on the dark couphng ^ presented in 
Ref. [16], adding to the analysis the contribution from the expansion rate pertur- 
bation /C and imposing adiabatic initial conditions for all fluids. We have therefore 
modified the Boltzmann CAMB code [36] to incorporate the dark coupling ^ and the 
IC terms. 

In the synchronous gauge, /C = 9T/{S'H)+h/{6'H) and the perturbation equations 
reduce to: 



■'dm 



•^ Pdm Pdrr 



kVT h 

-+6 



Vdr, 



-Hu 



dm 



Sde = -(1 + w){kVde + -h) - 3H (1 - W) [5de + V- (3(1 + w) + ^^^ 



Vde 



2H{1 + 



e 



l + w 



Vde + :r- Ode - Kri 



l + w 



l + w 



(4.1) 

(4.2) 
(4.3) 



(4.4) 



where vt is defined in Eq. ( [A.IOI) . 

We have extracted the cosmological parameters by means of the publicly avail- 
able Markov Chain Monte Carlo package cosmomc [37]. The cosmological model is 
described by ten free parameters 



{1^6, i^dm, OcMB, T, ^k, fu, Ul, ^, ''^s, ^s} 



*See Rcf. [18] for Qa = ±rpiim, where F is a constant. 
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where Ub = flbh'^ and Udm = ^dmh^ are the current baryon and dark matter densities 
respectively, Ocmb is proportional to the ratio of the sound horizon to the angular 
diameter distance, r is the reionization optical depth, VL^ is the spatial curvature, 
fy = Vty/VLdm refers to the neutrino fraction, n^ is the scalar spectral index and As 
the amphtude of the primordial spectrum. 

The analysis is restricted to negative couplings and also w > —1 (to ensure 
the avoidance of phantom behaviour), exactly as it we did previously in Ref. [16]. 
The basic data set we exploit here includes a prior on the Hubble parameter of 
72 ± 8 km/s/Mpc from the Hubble key project (HST) [38], the constraints coming 
from the latest compilation of supernovae (SN) [39], the matter power spectrum 
(large scale structure data or LSS data) from the spectroscopic survey of Luminous 
Red Galaxies from the Sloan Digital Sky Survey survey [40], the H{z) data from 
galaxy ages [41] and the WMAP7 data [1,35]. 

CMB constraints the amount of dark matter at redshift ~ 1000. In the presence 
of a negative dark coupling, the energy flows from dark matter to dark energy, thus 
dark matter energy density is smaller today as it can be seen in Fig. |^ (left panel). 
This effect is compensated for large scale structures by a larger growth of dark 
matter perturbation (see e.g. [42]). Figure |I|, left (right) panel illustrates the la and 
2(j marginalized contours obtained in the ^-Qdmh"^ {^^^) plane. We verifled that 
the results do not differ signiflcantly if including WMAP5 data (as we had done in 
Ref. [16]) instead of WMAP7 data. 

Overall, the results show that the addition to the analysis of the perturbation 
expansion rate K, leaves basically unaffected the quantitative constraints on the cos- 
mological parameters previously obtained in Ref. [16]. Indeed, all the additional 
terms introduced to make perturbations gauge invariant give negligible contributions 
at observable scales. 

5. Conclusions 

Interacting dark energy-dark matter cosmologies in which the coupling term is pro- 
portional to the Hubble expansion rate are revisited. While in previous works the 
perturbation in the Hubble expansion rate was neglected, it is illustrated here how 
the inclusion of such a term is mandatory to satisfy the gauge invariance of the the- 
ory. It also serves as a guide to deflne a covariant formulation of the dark sector 
interaction. In this work, the latter has been chosen to be expressed in terms of the 
expansion rate associated to the total fluid. This choice is however not unique, we 
could have used the expansion rate of any other fluid. For the case under study, we 
compute the linear perturbation evolution using a gauge invariant formalism. After 
imposing adiabatic initial conditions on the matter and radiation fluids, we flnd that 
the initial conditions for the coupled dark energy fluid are also adiabatic. This result 
is independent of the choice in the covariant formulation of the expansion rate. The 
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new terms arising from the expansion rate perturbation have neghgible quantitative 
impact on the constraints on cosmological parameters previously obtained in the 
hterature. A new analysis has been performed using the latest WMAP7 data. 
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A. Gauge invariant formalism 

The conventions we use are mostly from Ref. [22] with a few exceptions. For pertur- 
bations in flat space time, the perturbation variables can be expanded by harmonic 
functions Y^^\x,k) satisfying to (V^: + k'^)Y^^^ = 0. In the following we focus on 
scalar perturbations for which we define: 

if > = -iy^f , (A.i) 

A.I Metric perturbations 

For the metric defined in Eq. ( p.l|) , expanding in the Fourier basis the independent 
perturbations, we denote: 

A -^ ly(^) , 

H,^ ^ HLl^J + HtY^P , 

where Hij'-/^^ = 0. From now on, for sake of simplicity we will drop the tilde symbols. 
Remember that all these quantities are represented by the correspondent Fourier 
expansion and depend only on time and on the 3-momentum k, while the position 
dependence is left only in the basis Y elements. 

Gauge transformations are associated to infinitesimal coordinate transformations 
under which: [x^, x*) — )■ {x^, x*) = (x° — T, x* — L*). It can be shown that the metric 
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perturbation transforms as: 



B-B = -kT - L , 
Hl-Hl = Hl + kL/3 + UT 
Hx — Hx = Hx — kL . 



(A.3) 
(A.4) 
(A.5) 
(A.6) 



Before going to the gauge invariant variable definition, let us define some useful 
metric quantities and their transformations: 



(Tn 



^{HT-kB 

rC 



n = HL + -HT, 



/c 



H 



-HA + -vt + Hl 



where vt is the center of mass velocity for the total fluid, satisfying 

(1 + Wt)Vt = ^(1 + Wa)^aVa ■ 
a 

In the text is also sometimes used the following quantity: 



(A.7) 
(A.8) 
(A.9) 

(A.IO) 



/Ca 



1 
H 



-UA + -Va + Hl 



(A.ii: 



The physical meaning of the quantities above is the following: ag represents the 
shear perturbation, TZ is the curvature perturbation and /C (/Cq) is the expansion 
rate perturbation of the total (a) fluid. These quantities are not gauge invariant but 
transform as: 



^ff-S 



kT , 

n-n=n+HT 



}C-}C = ^('H-'H^]t = =T 

n\ J H 



(A.12) 
(A.13) 

(A. 14) 



where H = Ti/a is the usual Hubble parameter defined in the proper time. From the 
definition of Eq. ( |A.14| ) we see explicitly that we can identify /C as the perturbation 
oiH. 

We now define gauge invariant quantities associated to the metric and fluid 
perturbations. Bardeen metric gauge invariants are defined [24] as: 



^B = A-jag 



$ 



B 



1 

3' 



1 

n 



Hl + -Ht - j(Tg . 



(A.15) 
(A.16) 
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One can also build the following gauge invariant observable related to the expansion 
rate perturbation: 



It is also useful to define the following gauge-invariant quantity: 



A=<ilB-^-(l-^]^B = lil+WT){VT-^B 



(A. 17) 
(A.18) 



(A. 19) 



with Vt the (reduced) gauge invariant velocity of the total fluid defined by: 

{l + WT)VT = J2il + Wa)naVa. (A.20) 

a 

A. 2 Useful equations 

The perturbation equations for the metric can be derived from Einstein equations: 



$R + ^ 



B ^ ^B 






A;2 pj, 



fc2 



n 



n 



x^ 



(A.21) 



^B-^ = ~{l + WT)Vr = \Y.^l + w^)n^r^ iv=^\ , (A.22) 



3(1 + Wt) + |X2 Ea 3 (1 + Wa) ^a + f^^ 



(A.23) 



where we have defined x = k/Ti. From the previous equation one can obtain the 
following relation for the expansion rate perturbations: 



C 



y-2(l + -T) 



Vt 



Ca = yK-^(l + WT)V^T. 



(A.24) 
(A.25) 



For the sake of completeness we also provide the relation between the entropy per- 
turbation Ta, defined in Eq. (|2.7| ), and the sound speed in the rest frame of the fluid 
c|^ which is given by: 



w, 



Xa — [Csa '^Aa) 



A„ 



pa f^B Va 



Pa \'H k 



(A.26) 
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